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1 Introduction 

Let G = iy^ E) be an undirected graph without loops or multiple edges (called the host graph), 
with n vertices (its order), and m edges (its size). Our problem is stated as follows: 

Problem 1 (MaxDDBS) Given a connected undirected host graph G, an upper hound A for 
the maximum degree, and an upper hound D for the diameter, find the largest connected subgraph 
S with maximum degree < A and diameter < D. 

MaxDDBS is a natural generalization of the well-known Degree-Diameter Problem (DDP), 
which asks for the largest graph with given degree and diameter [8]. DDP can be seen as 
MaxDDBS when G is the complete graph K^ for sufficiently large n. Problem [T] was recently 
introduced in [1], where various practical applications are discussed and a heuristic approxima- 
tion algorithm to solve MaxDDBS is given, since it is computationally hard. 

Regarding computational complexity, MaxDDBS is known to be AAP— hard, since it contains 
other well-known AAP— hard problems as subproblems. In fact, restricting the search to only 
one constraint (either on the degree or the diameter), is enough to ensure AAP— hardness [4]. 
The Largest Degree-Bounded Subgraph Problem is AAP— hard as long as we insist that the 
subgraph be connected, but can be solved in polynomial time otherwise (Problem NDl of [3j). 
On the other hand, the Maximum Diameter-Bounded Subgraph becomes the Maximum Clique 
for D = 1, which was one of Karp's original 21 AAP— hard problems ^. MaxDDBS also 
turns out not to be in Apx, the class of AAP— hard optimization problems for which there is a 
polynomial-time algorithm with a constant approximation ratio. 

As mentioned above, MaxDDBS is also closely related to the Degree-Diameter Problem (DDP), 
stated by Elspas in 1964, which consists of finding the largest graph with a given maximum 
degree A and a given diameter D. Since the order of such a graph cannot exceed the quantity 

Ma,d = 1 + a + A(A - 1) H h A(A - 1)'^~\ called the Moore hound, if we take G as the 

complete graph on Ma,d vertices (denoted by Km^ b) ™ ProblemfTl we get the Degree-Diameter 
Problem. Note that this does not imply that DDP is AAP— hard; actually, the complexity of 



DDP is not known to-date. Obviously, the Moore bound is also a theoretical upper bound for 
MaxDDBS. 

A graph whose order is equal to the Moore bound is called a Moore graph. Moore graphs are 
very rare; they exist only for certain special cases; for diameter D = 1, Moore graphs are the 
complete graphs of order A + 1, for maximum degree A = 2, Moore graphs are the odd cycles. 
The only other Moore graphs are of diameter D = 2, A = 3,7 and possibly 57, [8j. We denote 
by N/\^£) the order of the largest graph that can be constructed with maximum degree A and 
diameter D; the current lower bounds for A^a.d are shown in [6J. 

A case of special interest is when the host graph G is a common parallel architecture, such as the 
mesh, the hypercube, the butterfly, or the cube-connected cycles. If there are any constraints 
on communication time between two arbitrary processors, then MaxDDBS corresponds to the 
largest subnetwork that can be allocated to perform the computation. The case of the mesh and 
the hypercube as host graphs were already treated in [1] , where some bounds were found for the 
order of MaxDDBS in a A:-dimensional mesh. Here we revisit in more detail the case of the 
mesh as a host graph. We refine the bounds given in [Ij for the order of the largest subgraph in 
arbitrary k > 1, and we focus on the cases A; = 3, A = 4 and A: = 2, A = 3. For those particular 
cases we give constructions that result in larger lower bounds. 

The rest of the paper is organized as follows: In Section [2] we give the bounds for MaxDDBS 
in a /c-dimensional mesh. Then, in Section [3] we give the constructions for A = 4 in the 3- 
dimensional mesh. Section |4] gives constructions for A = 3 in dimension two. Finally, in Section 
[5] we discuss some open problems and research directions. 

2 MaxDDBS in the /c-dimensional mesh 

Here we will assume that the host graph G is an infinite A;- dimensional mesh, and we are looking 
for a subgraph of maximum degree A < 2k, and diameter D. We can associate our mesh with 
an L^ metric space in dimension k. Pick an arbitrary point in this A;-dimensional L^ metric 
space as the center of a coordinate system. Now, the vertices of the mesh are the points with 



integer coordinates {lattice points), and two lattice points are joined iff they are at distance 1. 
The largest subgraph of degree A = 2k and diameter D corresponds to a closed ball of radius 
D/2. The number of lattice points contained in a ball of radius D/2 is variable, and depends on 
the location of the center of the ball. The maximum number of lattice points is achieved when 
the center of the ball is a lattice point itself, for even D, and when the center of the ball is the 
midpoint between two adjacent lattice points, for odd D. Such balls will be called maximal, and 
will be denoted Bk{p), where p = [^J i depending on its parity, and A; > 1. Figure nl depicts two 
maximal balls in dimension two, with diameters 5 and 6, respectively. 
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Figure 1: Maximal balls in the two-dimensional L^ metric space 



For more details about the shape and location of maximal balls, we refer to the comprehensive 
study by Dougherty and Faber [2, where the same problem appears in a slightly different 
context. In that paper, the infinite /c- dimensional mesh is interpreted as the Cayley graph of 
the free abelian group Z'^, and the points with integer coordinates contained in the ball Bj.{p) 
correspond to words of length at most D in the canonical generators of Z . The aim of [2j was 
to construct large Cayley graphs on abelian groups, with given degree and diameter. 

In order to simplify notation, we will also use Bk[p) to denote the set of points with integer 
coordinates contained in the closed ball Bk{p). The order of the largest subgraph S of degree 
A < 2k and diameter D = 2p or D = 2p + 1, that can be constructed on the A;-dimensional 
mesh, will be denoted Ni:{A,p). Alternatively we could use the notations Bk{D) and Ni^{A, D), 
specifying whether D is even or odd. If k' > k, the following inequalities are straightforward: 



\Bkip)\<Nki2k,p)<\Bk>ip)\ 



(1) 



The first inequality tells us that if we go to a higher dimension, keeping A and p constant, 
we can construct larger subgraphs. The reason for that is that we can move along the extra 
dimensions in order to avoid 'collisions'. Figure [2] is an example of one such construction in 
dimension A; = 3, of a subgraph with degree A = 4, diameter D = A (i.e. p = 2), and 18 vertices, 
whereas |i?2(2)| = 13. 



Figure 2: Construction for A = 4 and Z) = 4 in the 3D mesh 



From (II| follows the importance of determining the numbers |-Bjt(p)|. Counting the number 
of lattice points contained in a circle is a problem that goes back to Gauss, and there are 
several approximate results for the number of lattice points in balls and other sets, e.g. |121ll5j. 
Regarding the exact number of lattice points contained in closed balls in the L^ metric in 
arbitrary dimension, the main reference seems to be a paper by Vassilev-Missana and Atanassov 
The following result was given in p[], and we reproduce it here with slight modifications. 



Theorem 2.1 The cardinality of Bk{p) is 
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Proof: 



The ball Bk[p) can be constructed as the union of a smaller ball in the same dimension /c, plus 
two balls in dimension k — 1. Let D be even, and let us place the origin of our coordinate system 
in the central lattice point. Then the subset consisting of all the lattice points having the /c-th 
coordinate x^ equal to zero is B^^iip). This subset separates B).{p) into two hemispheres, one 
made up by those lattice points with a positive k-th. coordinate, and those with a negative fe-th 
coordinate. The layers with Xk = ±1 are Bk-i{p — 1), i.e. they have diameter D — 2 = 2{p — 1). 
If we remove one of these layers (say, the one with Xj. = —1), and put both hemispheres together, 
we get Bk{p — 1). Figure [S] shows the decomposition for A; = 2 and p = 3. The case of odd 
diameter is very similar. 



*- — « ^^ 

•^ 1( (I II ;• 

• " — ;*' 



Figure 3: Decomposition of the ball of diameter 6 



From this decomposition we get the recurrence relation 



f{k,p) = fik,p - 1) + f{k -l,p) + f{k - l,p - 1) 



(3) 



where f{k,p) denotes the number of lattice points in Bk{p). The boundary conditions are: 



f{K^) 



1 if D is even 

2 if L> is odd 



fi-^.P) = { 



(4) 



2p +1 if D is even 
2(p + l) if D is odd 

We want to find the generating function Ak{z) = ^^pX) f{k,p)zP. Multiplying ([3J) by z^ and 
summing over p > 1 we get 



Ak{z) - Ak{0) = zAkiz) + (Ak-iiz) - Afc-i(O)) + zAk-iiz) (5) 

whence 



Ak{z) = ^Afe_i(z) (6) 



With the aid of the boundary conditions Q we get 



jl+z)'' 



iiD = 2p 



M^)={ ^'~t;l ' (7) 

2ji±fL_ ifD = 2p+l 



For even D, Ak{z) is the product of (1 + z)'' = Y,p (p)-^^ and 1/(1 - zf+^ = J2p C~^^)zP- Then, 
the series of Ai.{z) can be obtained as the convolution of the respective factor series. The series 
of Ak{z) for odd D can be obtained in the same manner. D 

For even D, the numbers \Bk{p)\ turn out to be the Delannoy numbers (sequence A008288 of 
[9]), which appear in a variety of combinatorial and geometric problems |13]. This particular 
interpretation of Delannoy numbers was first given by Vassilev-Missana and Atanassov [14J , and 
later rediscovered by Schroder |llj . and then by us. Our formulation and proof are different 
from the ones in |14^ lllj . For odd D, the numbers \Bk{p)\ are known as a Riordan array of 
coordination sequences (sequence A113413 of [9]). Tables [l] and ^ show the first few values 
of \Bf^{p)\ for even and odd D. They can be constructed in a Pascal-like fashion, with the 
convention that |i?o(p)| = 1- 

It is known that Delannoy numbers have no closed form, meaning that they cannot be represented 
as a linear combination of a fixed number of hypergeometric terms (which can be verified with 
the aid of the methods developed in [lO]). However, we can extract asymptotic information from 



the generating function Ak{z) in the proof of Theorem 2.1 above. Recall that a is an algebraic 



singularity of the function f if f can be written near a as 



^(^) = A(-) + (T^ W 
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Table 1: Some values of \Bk{p)\ for even D 
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Table 2: Some values of |i?fc(p)| for odd D 



where /o and g are analytic near a, g is nonzero near a, and w is a real number different from 
0,-1, —2, .... We readily recognize that 1 is an algebraic singularity of j4fc(z), since Af^iz) can 
be written in the above form, with g{z) = (1 + z) for D even, and g{z) = 2(1 + z) for D 
odd, and all the other conditions are satisfied. Now we can readily apply Theorem 3 of [7J: 

Theorem 2.2 Suppose that for some real r > 0, A{z) is analytic in the region \z\ < r, and has 
a finite number t > of singularities on the circle \z\ = r, all of which are algebraic. Let Oj, Wj, 
and gi be the values of a, u, and g in IM), corresponding to the i-th such singularity. Then A{z) 
is the generating function for a sequence {an) satisfying 

where Vt is the maximum of the uoi and T denotes the Gamma function. 
We get 
Corollary 2.1 

3 Subgraphs of degree 4 in the S-dimensional mesh 

An interesting special case is /c = 3 and A = 4. In this case, the inequalities M translate to 
\B2{p)\ < N^{A,p) <\B^{p)l OT 

2p2 + 2p + 1 < iV3(4, p) < 4p3/3 + 2p2 + 8p/3 + 1 ii D = 2p 

2(p2 + 2p + 1) < iV3(4, p) < Ap^/3 + 4p2 + l4p/3 + 2 ii D = 2p + I 

The following theorem shows that the lower bounds are in fact a lot closer to the upper bounds: 



Theorem 3.1 

V/3 + 2/ - 4p/3 + 3 < iV3(4,p) < 4p3/3 + 2p2 + 8p/3 + 1 
V/3 + V + 2p/3 < N3{4,p) < V/3 + V + l4p/3 + 2 

Proof: 



ifD = 2p 
ifD = 2p + l 



(9) 



Let D = 2p, with p > 2, and let us go back to our L^ metric space model in dimension two. 
W.l.o.g. we can pick the center of our coordinate system as the center of all our balls and 
constructions. Now let E2{p) be the graph obtained from i?2(p) by removing all the edges along 
the y-axis, and the two vertices that are left isolated (see Figure H|. 
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(c) 



(d) 



Figure 4: Construction for A = 4 and even diameter in the 3D mesh 

Note that in E2{p) there are two classes of nodes: the ones along the y-axis, shown as black 
squares, and all the other ones, shown as black circles. The interior circular nodes have degree 
4, while the interior square nodes only have degree 2, and we will use them to move along the 
third dimension z. 

Our three-dimensional graphs H3{p) (with p > 2) will consist of layers of £2(1) connected via 
the square nodes: At 2: = we have E2{p), and at z = ±i we have E2{p — i) , ior 1 < i < p. It is 
straightforward to check that H3{p) has diameter D = 2p, so let us concentrate on the number 
of nodes. 

The number of vertices on each layer E2{i) is |i?2(^)| minus two 'missing' vertices. Therefore 
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\E2{p)\ 
\H3ip)\ 



2p^ + 2p-l 



2 E?=i^(2^^ + 2i - 1) + (2p2 + 2p-l) 
= 4p3/3 + 2p2 _ 4p/3 + 3 

Now let D = 2p + 1, with p > 2. The construction here is also made with layers that are i?2(p) 
with the central 'spine' suppressed (i.e. the edges on the y-axis and the two tip vertices). Let us 
call these graphs 02{p)- Figure [s] shows 02(1)1 02(2), and 02(8). The three-dimensional graph 
Qz{p) has a layer of type 02(p — i) at z = iti, for < i < p — 1. 



(a) 



II ■ II 

II ■ 1) 



(b) 



(c) 
Figure 5: Construction for A = 4 and odd diameter in the 3D mesh 

Again, it is easy to check that the diameter is correct. Regarding the number of nodes, in every 
02(2) there are only two nodes missing from B2{i). Hence 



\QM\ =4E£i(i' + 2i) + 2p2 + 4p 
= (4p3 + 12p2 + 2p)/3 



D 



Note that the above constructions are asymptotically optimal, since they agree up to the second 
term with the upper bounds. They are also optimal in another sense, as shown by the following 



Corollary 3.1 The average degree S of the constructions in Theorem 3.1 tends to 4 as p ap- 
proaches infinity. 
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Proof: 

6 = 2A/V, where A is the number of edges and V is the number of vertices. The number of 
edges is 8p^/3 + 4p/3 + 2 for I) = 2p, and 8p^/3 + V - 2p/3 for D = 2p+l. D 

It is also very hkely that these constructions can be extended to higher dimensions. 

4 Subgraphs of degree 3 in the 2-diniensional mesh 

The smallest case oi A < 2k that makes sense is A = 3 in dimension two. In this case we will 
also show that the lower bounds are quite close to the upper bounds. We have the following 

Theorem 4.1 

2p2 - 2p + 1 < N2{3,p) < 2p2 + 2p + 1 ifD = 2p 

8r2 + 2r = 2p2+p < N2{3,p) < 2p^ + Ap + 2 if D = 4r + 1 (10) 

8r2 + 10r + 6 = 2p2+p + 3 < N2{3,p) < 2p'^ + Ap + 2 if D = Ar + 3 

Proof: 

As in Section [3j we will give two constructions that achieve the new bounds. In this case, a 
rigorous mathematical description of our constructions will be quite cumbersome; instead, we 
will use some geometric analogies to describe them. What we do is that we take the balls B2{p) 
and we try to fill them up in a convenient way with as many 'building blocks' as possible. Our 
construction elements will be 2 x 2 square blocks, and 2x1 half-blocks, or rectangular 'bricks'. 
We have to pack them up in such a way that no four edges meet at a single point. 

Let us start with even D: We place as many square blocks as possible along the x and y axes, 
and if there is soine space left at the end, we fill it up with bricks, as shown in Figure ro] (a). 
On each horizontal semi-axis we will have to use {p — l)/2 square blocks, and on each vertical 
semi-axis we will need {p — 3)/2 square blocksrl 

Now we will be left with four empty triangular regions that will have to be occupied. They are 
isosceles triangles, with side p — 2. We will fill them with bricks, forming a pattern of interwoven 



^A fraction here means that we have to complete with a brick at the end. 
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Figure 6: Constructions for A = 3 in the two-dimensional mesh 

horizontal and vertical bricks, as shown in Figure [7j Starting from the innermost corner there 
is only one way to do this, since every other vertex on the triangle sides already has three edges 
attached to it. Moreover, all triangles constructed this way will be the same, save rotations and 
reflections. We will make use of this fact again in the case of odd D. The construction ensures 
that no vertex has degree 4. 




Figure 7: Brick pattern in a triangular region 

Next we have to show that the diameter is D = 2p. First consider the blocks along the axes. We 
shall call the vertices belonging to one of those blocks 'block vertices'. It is clear that starting 
from any block vertex we can reach any other block vertex in at most 2p steps. Hence we only 
have to include the vertices in the triangular regions, and show that for any such vertex, we can 
reach any other vertex in at most 2p steps. The procedure is quite simple: First note that we 
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can reach the innermost vertex of the triangle (which is also a block vertex) in at most p — 2 
steps (the length of the triangle sides). Call this vertex a 'corner'. From there we can reach any 
block vertex in at most p + 2 steps, and we can reach any other corner in at most 4 steps, which 
means that we can reach any other vertex located in a triangular region in at most 2p steps, as 
desired. 

Now we are going to count the number of vertices: If p is even we have 2(p — 3) square blocks, 
plus four bricks, one at the end of each semi-axis. The square blocks contribute 10{p — 3) + 3 
points, and we have to add 12 additional points corresponding to the four bricks. Each triangle 
contributes 2[{p - 2)/2j2 - 2[{p - 2)/2j = [{p - 3f - l]/2 points. Hence, the total number of 
vertices is 2p^ — 2p-\-l. li p is odd we have 2{p — 2) square blocks, and no bricks at the end of the 
semi-axes. The square blocks contribute 10{p — 2) -|- 3 points. Finally, each triangle contributes 
2[(p — 2)/2j^ = {p — 3)^/2 points. Hence, the total number of vertices is also 2p'^ — 2p -|- 1. 

Now consider the construction for odd D. In this case we have a central 'spine' of length 
D = 2p + 1 (the y-axis), and two horizontal semi-axes of length p. We place \{p — 2)/2] square 
blocks on each horizontal semi-axis, as shown in Figure M (b). Again, we are left with four 
triangular regions that are filled up with interwoven bricks, in the same manner as before. In 
this case, two triangles have sides of length p, and the other two have sides of length p — 1. 

We show that the diameter is -D = 2p -|- 1. We proceed as in the case of even D. In this case 
we distinguish block vertices and 'axis vertices' (i.e. vertices belonging to the central spine or 
to the horizontal semi-axes). We also distinguish two central vertices at the intersection of the 
vertical spine and the two horizontal semi-axes. It is clear that starting from any block or axis 
vertex we can reach the closest central vertex in at most p steps. Now, if we start at a vertex 
belonging to one of the larger triangles, of side length p, we can reach the corresponding corner 
(which is one of the central vertices) in at most p steps. Starting at a vertex in one of the smaller 
triangles, of side length p—1, we can reach the corresponding corner in at most p — 1 steps, and 
the nearest central vertex is just one step away. Since both central vertices are adjacent, the 
proof of the diameter follows. 
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Regarding the number of vertices, we have 2p + 2 vertices on the vertical spine, plus 2p vertices 
from the horizontal semi-axes, plus 6\{p — 2)/2] block vertices, plus the triangles. If p is even, 
the two large triangles contribute p{p — 2) extra vertices, and the two smaller triangles contribute 
{jp — 2)^ extra vertices. If p is odd, the two large triangles contribute {p — 1)^ extra vertices, 
and the two smaller triangles contribute {p — \){p — 3) extra vertices. Hence the total number 
of vertices is 2p'^ + p \i p is even, and 2p'^ + p + 3 if p is odd. 

D 
These constructions are also asymptotically optimal in the sense of the average degree: 



Corollary 4.1 The average degree 5 of the constructions in Theorem 4-1 tends to 3 as p ap- 
proaches infinity. 

Proof: 

The number of edges is 

3p^ -Qp + 6 a D = 4r 

3p^ -p-1 a D = 4r + l 

(11) 
3p2 - 6p + 4 a D = 4r + 2 

3p^ -p + 3 a D = 4r + 3 

D 

Despite being asymptotically optimal, these constructions can be improved without much effort 
for particular values of D, but extracting a general pattern and counting the vertices for arbitrary 
D may prove a challenge. In Table [3] we give the order of the largest graphs that have been 
constructed for some small values of D. The actual graphs can be seen in [6]. Note that for 
diameters D = 2, . . . , 6, the upper bound given in the table is smaller than the upper bound 



given by Theorem 4.1 and the largest graphs are optimal. These graphs have been obtained by 



exhaustive computer search. 
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Diameter 


2 3 4 


5 


6 


7 


8 


9 


10 


11 


12 13 


14 


15 


16 


Upper bound 
Order of largest graph 


4 6 10 
4 6 10 


14 
14 


22 
22 


32 

28 


41 
37 


50 

44 


61 
52 


72 
68 


85 98 
77 90 


113 

104 


128 
124 


145 
135 



Table 3: Orders of the largest known subgraphs of the two-dimensional mesh, with A = 3 and 
2 < -D < 16. 

5 Conclusions and open problems 

To us, it is always amazing that a simple combinatorial setting like this one provides such a wealth 
of interesting and difficult problems. In the preceding two sections we have seen constructions 
that result in lower bounds for Nk(A,p), that are asymptotically close to the upper bounds, for 
small values of k and A, with A < 2k. Now, is that also possible for arbitrary k and A < 2k? 
How close can we get to the upper bounds? The constructions in Section [3] could be generalized 
to higher dimensions, more precisely to the case A = 2A; — 2. The constructions in Section [4] 
might also be extended to A = 5 in dimension 3, but it does not appear to be easy. 

There are several additional properties that could be considered for every construction, like 
connectivity {sic fault-tolerance), average path length, symmetry, etc. For example, in parallel 
computing applications, our constructions of Section |3] might not be the best ones, as they 
impose a relatively high communication overhead on the square vertices. 

On the other hand, what is the computational complexity of finding the largest degree&diameter 
bounded subgraph in the mesh? Does it remain AAP— hard for all dimensions?, for some dimen- 
sions?, for any dimension? If it remains AAP— hard, then, can it be approximated to within a 
constant ratio? 

This same study can be carried out for other host networks of theoretical and/or practical 
importance: the hypercube, the butterfly, the cube-connected cycles, Cayley graphs, etc. In 
[1] there is a preliminary discussion for hypercubes, and a very incomplete heuristic study was 
attempted for some random networks, but other than that, the field remains totally virgin. 

Finally, it is worth noting that this problem could be related to a well-known problem arising in 
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parallel computing, namely that 'embedding' or 'emulating' an arbitrary graph in a mesh (see 
|16j . for instance). This connection was indicated to us by one of the reviewers, and it is also 
interesting to explore further. 
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